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dependence on period 
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critical current density measurements also show a scaling that indicates the super- 
conductivity properties are not dependent on the number of bilayers. 
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I. INTRODUCTION 



The level of interest in the superconducting proximity effect for layered systems has 
continued to be very high with the advent of structures using ferromagnetic materials in 
proximity with superconductors directly or via normal metals^""-. Theoretical studies of 
proximity systems continue to be carried out as predictions of the superconducting transi- 
tion temperature of bilayers and trilayers are critical in their use as detectors, as well to 
understand the interactions of the layers.-"- 

The dependence of the critical temperature as the number of periods in a proximity system 
is varied has also been recently considered.™ In this work, the system N[SN]™ was studied in 
the Nb/Cu system, showing how the transition temperature varied as the number of periods, 
m, were increased. Here N refers to a normal metal layer and S to a superconducting layer. 
An open question that still remains in light of this extensive work is exactly how does the 
critical temperature vary in the system [NS]™"? The open question of how a system would 
transition from a single bilayer to a multilayer is still not fully understood, as discussed under 
the DeGennes-Werthamer (DGW) model,- which is equivalent to the single mode Usadel 
model. In the work by Broussard,-ii^ the implication is that the transition temperature of 
a [NS]™ system would not change as m is varied, even though it would change for both a 
N[SN]™' as well as a S[NS]™' system. In a study of the Nb/Pd multilayer system, this result 
was used to calculate the transition temperature successfully, but no dependence on m was 
studied.— 

In an attempt to bring some clarity to this question, we have carried out a study of the 
Mo/Nb layered system, growing the [NS]"^ system with m varying from 1 to 4. The Mo/Nb 
multilayered system has not been previously studied, although Mo/V has been studied 
extensively by Karkut et alM In that work, anomalies in the superconducting transition 
temperature were observed and partially linked to the properties of the vanadium layers. 
In this work we have chosen the Mo/Nb system partly due to it not being studied before 
and the ability to see if what was observed for the Nb/Zr system^ would be replicated in 
a different system. In addition we have carried out measurements of the inductive critical 
current (current in plane) as a function of temperature to see how this varies as the number 
of periods is changed, since as pointed out in the work on critical currents in Nb/Fe bilayers 
by Geers et al.^^ that while transition temperatures probe the maximum value of the order 
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parameter, critical currents involve an averaging over the entire layer thickness. 



II. THEORETICAL CONSIDERATIONS 

Initially, it was believed that increasing the number of periods in a [NS]™ system would 
eventually be equivalent to an infinite multilayer,— which would imply a reduction in Tc as m 
increases, since an infinite multilayer is equivalent to a bilayer with layer thicknesses reduced 
by 1/2. Other work- however implies that the superconducting transition temperature (T^) 
for the [NS]™ layered system should not change as m increases. While it may initially seem 
counterintuitive that boundaries would be significant when approaching an infinite number 
of layers, viewing the system as analogous to a one-dimensional quantum potential well 
does lead us to expect this behavior.J^ In this analogy, each nonsuperconducting layer N is 
treated as a potential barrier and each superconducting layer S as a zero potential region. 
The amphtude of the wavefunction is related to the density of the superconducting pairs in 
the multilayer. 

To compare to our measurements, we use the dirty limit DeGennes-Werthamer theory 
used in previous studies of bilayers.— We use the formulation as laid out in previous work.- 
Calculations can easily be carried out for our system, with all parameters measured except 
for the electronic specific heat coefficients. For these we use published values of 7.8 and 2.0 
mJ/(mole-K^) for Nb and Mo, respectively.— 

As mentioned earlier, there has been recent work^ on how the Tc will vary for a N[SN]™' 
using the Usadel theory. For this configuration of layers, the Usadel model predicts that 
Tc increases as m increases, which is what was also seen in the earlier modeling under 
DeGennes-Werthamer theory for a N[SN]™ system.- Since the Usadel model in the single 
mode limit and the DGW theory are the same, it would be important to know for our films if 
a single mode limit is appropriate. Comparing the simple DGW model to the more complex 
Usadel predictions for the data of Kushnir et al.^ in Fig. [T]we find that the two agree within 
0.2K for the niobium thickness range we are using. In general the two sets agree over a wide 
range of niobium thickness. The lack of exact agreement at the extremes of the plots are due 
mainly to the DGW model having fewer fit parameters. For the calculation shown in Fig. [H 
there are no adjustable terms, as all have come either from the measured values in the paper 
by Kushnir et al. or in the case of the electronic specific heats, from Kittel.— Considering 
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FIG. 1. Comparison between a DeGennes-Werthamer model and the Usadel model for the 
Cu/Nb/Cu trilayer set CSl in Kushnir et alJ . For the DGW model, the resistivity of the lay- 
ers was taken from Kushnir et al. and the electronic specific heat coefficients were taken from 
Kittel.-il 



this limitation, the agreement is encouraging. Thus we feel that using the DGW model here, 
as it can be extended to multiple layers, is justified. The Usadel model, as developed for 
multiple periods in a N[SN]'" system depends on a center of symmetry,- which the [NS]™ 
system does not have. 
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III. SAMPLE PREPARATION 



Our multilayers were prepared using DC magnetron sputtering in a high vacuum system. 
The background pressure, P^, of the sputtering chamber was consistently less than 27/iPa 
before deposition. Silicon substrates were clamped to a room temperature platform, which 
rotated at 10 rpm to ensure a uniform sample deposition distribution onto the substrates. 
Ultra high purity (99.999%) Argon was used as the sputtering gas, with a pressure of 0.27 
Pa during sputtering. Separate DC power supplies operating at 250 W constant power 
were used for the Nb (S) and Mo (N) sources, with each sputter target of material rated at 
99.95% purity. Each source had a computer controlled mechanical shutter which could be 
closed to prevent the source from being deposited. The deposition rates for Mo and Nb were 
approximately 0.22 and 0.16 nm/s, respectively. Using these rates, we grew 36.9 nm Mo and 
43.2 nm Nb layers. Samples with 1-4 repeated bilayers were grown in different depositions. 



IV. X-RAY CHARACTERIZATION 



In order to characterize the lattice structure and spacing of our samples X-ray diffraction 

data along the sample normal was collected, for individual layers and the two period bilayers 

with an example shown in Fig. |2l The samples appear to have (110) texture, as is expected 

for samples of Nb and Mo grown at room temperature. Gaussian curves were fit to the 

diffraction peaks to give the Bragg angle and the width of the diffraction peak. The planar 

spacing (d) and Debye-Scherrer length (t) were then extracted according to the equations 

^_ A 
2sin{9B) 

and 

t - 

cos{9b)w ' 

where A is the wavelength of the X-ray beam, equal to 0.1542 nm, 6b is the Bragg angle, 
and w is the full width at half maximum. The rf(iio) was 0.2343 nm for Nb and 0.2229 nm 
for Mo, which points to the niobium layer being slightly stressed, while the molybdenum 
layer is nearly at the bulk value. The Debye-Scherer lengths are approximately 17 nm and 
24 nm for niobium and molybdenum, respectively, both less than the layer thickness, but 
we notice that as the molybdenum layer is under less stress, the correlation length is also 
longer. We see no evidence of multilayer peaks at these large values of layer thicknesses. 
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FIG. 2. X-ray diffraction along the sample normal for a two period (Mo/Nb) layered sample 
showing the (110) niobium and molybdenum lines. 

V. TRANSPORT MEASUREMENTS AND SUPERCONDUCTING 
TRANSITION TEMPERATURE 

To measure the transition temperature of each of the films, we first measured the electri- 
cal resistivity of the samples using the van der Pauw technique.— Contacts were made with 
pressed indium on the perimeter of the sample, which was then pressed into an indium foil 
sheet against the OFHC (Oxygen Free High Conductivity) Cu sample stage. Cooling was 
provided by a closed cycle cryocooler capable of reaching 6 K. The resistance data was gath- 
ered as the sample cooled from room temperature. Once the material went superconducting, 
we warmed and cooled it around the transition temperature to establish consistency. These 
measurements at low temperature were repeated for both of the van der Pauw configura- 
tions. The temperature at which the transition occurred is taken to be the midpoint of the 
transition; the width of the transition (from 10% to 90% of the normal state resistivity) is 
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typically O.OIK. 

Pure films of Nb and Mo of approximately 90 nm were grown and showed low residual 
resistivity ratios (RRR), indicating the samples are in the dirty limit of superconductivity. 
For Nb the RRR was 3.0 and low temperature resistivity(pioK) was 69 nfi-m, while for 
the Mo film RRR was about 1.7 and pioK was 56 nfi-m. The of our pure Nb film was 
8.77 K. We were unable to tell if the Mo film was superconducting, so we will consider it 
to be a normal metal for our purposes. (It is known that the of Mo is very sensitive to 
contamination and easily suppressed below the bulk value of 0.9 K.) 

The electrical resistivity for the Mo/Nb layered films showed a RRR between that of Nb 
and Mo, at a value of 2.7 for the 4 period sample decreasing to 2.3 for the one period sample. 
In an ideal situation, we would not expect RRR to change, but we did see both it and the 
thermal component of the resistivity (p290K — Piok) changing, the latter increasing as the 
number of periods decreased. The low temperature resistivity of the layered structures was 
of order 50-80 nfi-m. 

To verify the quality of the samples, we also measured the superconducting transition 
using an inductive coil technique.^- A coil with 64 turns of 1/4 mm diameter wire with coil 
inner diameter 3.0 mm and 3.5 mm outer diameter was pressed against the film surface, and 
excited with a 100 kHz drive current of 0.1 mA rms. The out of phase signal was measured 
using a lock-in method. In Fig. [3] we show the comparison of the transitions for the resistive 
and inductive measurements and see as expected that the inductive transition starts as the 
resistive transition is ending. The inductive transitions are also typically 0.01 K wide. 

VI. DEPENDENCE ON PERIOD 

As stated before,- what has been expected for a [NS]"* system as the number of periods 
(m) gets larger is for the system to go from the single to infinite case, while the theoret- 
ical investigation implies that the DeGennes-Werthamer theory would give an unchanging 
transition temperature as the number of periods increases. In Fig. H] we show the resistive 
transition temperature for the Mo/Nb layered sequences, where the number of bilayers goes 
from 1 to 4. Instead of the transition temperature decreasing as m increases, we see what 
looks like an increase going from m=l to m =2 with the values holding constant from then 
on. The Tc for the m=l agrees very well with the prediction from the DGW model. The 
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FIG. 3. Comparison of the superconducting transition for a 3 period Mo/Nb sample. The resistive 
transition is shown by the fihed circles and the inductive transition is shown by the filled squares. 

value for the infinite multilayer result would be about 1 K lower, which is off the graph. 
Clearly we do not see any decrease in the transition temperature towards the infinite multi- 
layer result, which is very similar to what was observed in 1-3 period Nb/Zr layered samples^ 
and continues to match those predictions. The increase seen in this figure for m > 1 is not 
similar to the work by Kushnir et ai, who predicted a significant increase in Tc (~ 0.5 K) 
as m increased from 1 to 2. In that work, since the fraction of superconductor is increasing 
with m, there is a large and sustained increase in T^, as m increases. We only see the jump 
up at m = 2. This could simply be explained by the niobium layers improving their Tc 
(due to both surfaces being protected from oxidation) by merely 0.15 K. Since the transition 
temperatures remain nearly constant for m > 2, we feel this is more in agreement with the 
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FIG. 4. Measured resistive superconducting transition temperatures for (Mo/Nb)m sequences vs. 
m the number of bilayers. In addition, the theoretical prediction for a bilayer is displayed. The 
value for an infinite multilayer is much lower, at 6.91 K for this system. 



[NS]*" system.-^ 



VII. INDUCTIVE CRITICAL CURRENTS 

To further probe the superconductivity of the multiple period structures, we measured 
the inductive critical currents in the films.^iS We use a 300 turn coil (made of 25 fim wire 
with inner and outer radial dimension of 1.05 mm and 1.95 mm, and an axial height of 0.46 
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FIG. 5. 3rd harmonic voltage vs. coil current for a pure Nb film at 7.38 K, showing how the critical 
current is defined 

mm) pressed against the samples, driven at 1 kHz and look for the nonlinear response at 
3 kHz. This obviously gives a current in plane type measurement. The coil specifics and 
definition of the critical current density are identical to that in Classeen et ai, as shown 
in Fig. O We point out that inductive critical current density measurements do not suffer 
from some of the issues that transport measurements do such as the pile up of current on 
the edges of the patterned film.— 

The measured critical current density as a function of temperature for the samples, shown 
in Fig. El follow a standard Ginzburg-Landau dependence, Jc cx (1 — t)^/^ where t is the 
reduced temperature, as expected. Our pure niobium sample has a value of Jc(0) that is an 
order of magnitude lower than predicted by the GL depairing formula in Geers' paper. This 
result might be due to vortex formation, since the induced current path in the film will be of 
the order of the width of the coil, which is much larger than the coherence length. What is far 
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FIG. 6. Log-log plot of the inductive critical current density vs {1 — t) where t is the reduced 
temperature for the pure niobium sample and the layered samples. The line gives the expected 
Ginzburg-Landau dependence. 



more relevant to our argument here is that the while the layered samples have lower critical 
current density compared to the pure niobium sample, the critical current densities for the 
layered samples collapse to essentially the same curve. As stated in Geers, since critical 
current density reflects an average over the layers, then this average is also unchanging as 
the number of periods varies. This gives even greater weight to the point that the order 
parameter shape is somehow remaining, on average, unchanged as the number of periods 
increases, which again goes against the usual perspective on how one would interpret the 
transition to an infinite multilayer result. 
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VIII. CONCLUSIONS 



As has been seen before, the superconducting transition temperature for a repeating 
bilayer structure composed of normal and superconducting layers seems to remain constant 
as the number of periods goes from 1-4, which is in agreement with previous theoretical 
analysis but is non- intuitive. The inductive critical current density is also seen to not depend 
on the number of periods. Obviously the behavior for larger m remains to be studied as 
well as other geometries, but this indicates work still needs to be done to fully understand 
the behavior of such systems with free boundary conditions. 
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